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We present a theoretical description of dissociative recombination of triatomic molecular ions 
having large permanent dipole moments. The study has been partly motivated by a discrepancy 
between experimental and theoretical cross-sections for dissociative recombination of the HCO"'" 
ion. The HCO'*' ion has a considerable permanent dipole moment {D «4 debye), which has not 
been taken explicitly into account in previous theoretical studies. In the present study, we include 
explicitly the effect of the permanent electric dipole on the dynamics of the incident electron using 
the generalized quantum defect theory, and we present the resulting cross section obtained. This 
is the first application of generalized quantum defect theory to the dissociative recombination of 
molecular ions. 
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I. INTRODUCTION 

Dissociative recombination (DR) of small molecular 
ions such as HCO+ or Hj^ plays an important role in the 
chemistry of interstellar clouds and is therefore one of the 
key elements in modeling the composition and temporal 
evolution of the clouds. In recent years, the theoretical 
description of the DR of Hj^ in collisions with low energy 
electrons [ll, 0, 0, 0] has exhibited good overall agreement 
with experimental data. Following the success of the the- 
oretical method developed for Hj^, a similar treatment 
was applied to describe DR in HCO+, motivated by 
the larger goal to develop a complete set of tools to treat 
dissociative recombination of small polyatomic molecular 
ions. 

The description of DR in HCO+ has proved to be far 
from simple. Although the recent theoretical studies 
0] have demonstrated that the non-Born-Oppenheimer 
Renner-Teller coupling plays an important role in DR of 
HCO+ (similar to the role of Jahn- Teller coupling in Hj^ 
DR), the theoretical cross-section obtained still remains 
about a factor of 2-3 smaller than the lowest experimen- 
tal result. Here, we point out that different experiments 
with IICO+ give quite different results for cross-sections 
and/or thermally-averaged rate coefficients of the process 
(see, for example. Fig. 3 of Ref. Q). Although the ex- 
perimental data differ from each other by up to a factor of 
4, it seems that existing DR theory is still missing some 
important effect(s). One of the possible improvements 
would be to account for the relatively large permanent 
electric dipole moment of 1100"*". The permanent dipole 
moment of HCO+ has been estimated in several studies 
0, H, [1] with values ranging from 3.3 to 4.0 debye (f de- 
bye is 0.393430 atomic units). In this study we adopt 
the value D = 3.9 debye from Ref. ^8J. Therefore this 
study is devoted to (1) the development of a theoreti- 
cal approach that accounts for the electron interaction 



with the electric dipole moment in addition to the usual 
Coulomb interaction between the ion and the electron 
and (2) an application of the approach to DR of IICO+. 
We will use atomic units (a.u.) in this study. 

One of the main theoretical techniques used in our 
previous studies is multichannel quantum defect theory 
(MQDT) [l3,[iil,[i3,[il- The standard version of MQDT 
[l0|,ll3| was initially developed for purely Coulomb poten- 
tials with centrifugal terms, i.e. potentials with asymp- 
totic behavior V{r) = — 7 + ^^77^, where r is the ra- 
dial electronic coordinates and I is an integer. A first 
generalization to the problem of a combined Coulomb 
and an attractive or repulsive dipole potential was in- 
troduced by Bely [3| in the context of electron scatter- 
ing from He"*", where the ionic dipole moment produces 
electron escape channels characterized by radial po- 
tentials having noninteger or even complex values of /. 
Refs. [11,13 

have generalized MQDT to account for the 
interaction between colliding partners behaving asymp- 
totically as with arbitrary a and positive p = 0,1, 
or 2. Subsequent studies [o, IH [H, [13, [U, IH have 
developed applications to a broader class of long-range 
potentials, e.g. with p = 3, 4 and 6 as well as to 
more general numerically-specified potentials. In par- 
ticular, the generalized theory (GMQDT) is well-suited 
for the present problem, where we want to account for 
the dipole-electron interaction in addition to the regu- 
lar Coulomb potential at large electronic distances. The 
GMQDT itself is developed in detail in Refs. [H, [H], 
and thus we only summarize here the main insights into 
the physics involved in main formulas of GMQDT, as 
well as describing the adaptations to the theory that are 
needed in order to describe DR of molecular ions with a 
permanent electric dipole moment. As is shown below, 
the new theoretical DR method to account for the dipole 
moment is very similar to the one developed earlier [1, [^ 
with comparatively small changes in the formulas for the 
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reaction matrix K and in the closed-channel elimination 
procedure. Therefore, a DR approach similar to the one 
developed in our earlier study Ref. Q is employed here. 
For this reason, in this communication we discuss only 
the main steps of the approach and the differences with 
the previous study 6]. 

Electron scattering from neutral polar molecules was 
treated in a frame transformation methodology by Clark 
[2^ and tested quantitatively by Clark and Siegel [13] ■ 
That theoretical formulation and the subsequent work 
of Fabrikant [25j bears similarities to the present one, in 
that the combined dipolar plus centrifugal potentials for 
the scattering electron are first diagonalized to find the 
scattering eigenchannels. 



II. HAMILTONIAN OF THE PROBLEM 

In order to model the DR process of HCO+, we write 
explicitly the Hamiltonian of the ion-electron system as 
H ~ Hicm + Hei, where Hion is the ionic Hamiltonian 
and Hei describes the electron and its interaction with 
the ion. Here, Hion is basically the vibrational Hamil- 
tonian of the molecular ion expressed in Jacobian coor- 
dinates Q={i?cOj -Rgh, <y5}, where i?co and Ran rep- 
resent respectively the distances C-0 and G-H (G is the 
center of mass of the CO bond), 6 is the angle between 
vectors OC and GH, ip is the angle of azimuthal ori- 
entation of the bending. Here, we consider Rgh as an 
adiabatic coordinate representing the dissociation path. 
The Hamiltonian Hion is explicitly written in Eq. (1) of 
Ref. In this approach, we neglected the rotational 
motion of the CO bond in space, but included relative 
rotation of H and CO. This approximation is justified 
by a large CO/H mass ratio. However, the approxima- 
tion may lead to an underestimated DR cross-section if 
the rotational autoionizing resonances in HCO play an 
important role in DR. As a result of the approximation, 
the projection rriip of the angular momentum L on the 
CO axis is conserved. We solve the Schrodinger equa- 
tion for the ionic Hamiltonian keeping the Rgh coordi- 
nate fixed: This determines vibrational wave functions 
^m^,i{RG}i',Rco,(^,'p) and corresponding adiabatic en- 
ergies Um^jiRoii) that depend parametrically on i?GH- 
Index I here numerates the ionic vibrational states for a 
given m,^. 

For the electronic part H^i of the total Hamiltonian H^i 
it was shown [2^ that only electronic states sa, ptt~ , per, 
and p7r+ have a significant contribution to the HCO+ DR 
cross-section. We include them into the treatment and, in 
addition, we also include the da states having a relatively 
small effect of DR. Therefore, the electronic Hamiltonian 
Hei ( Q) for a fixed value of Q and integrated over all elec- 
tronic coordinates can be written as an infinite number 
(Rydberg series) of matrices 5x5 (see Eq. (1) of Ref.jl]). 
In the matrix, the Renner- Teller couplings between pTT~ , 
pa and pn^ states are explicitly taken into account via 
the coupling coefhcients 7 and S. On the other hand, in 



our previous studies [1, Q , no coupling between sa, pa 
and da states was accounted for (except at short range 
in the ab initio calculation of the potential surfaces) be- 
cause these states are not coupled by the Coulomb field 
at large distances r in the basis of electronic states with 
a definite angular momentum. Now, if we include the 
electron-dipole interaction as well, it will introduce a 
coupling between states with different electronic angu- 
lar momenta. The coupling element between the sa, pa 
and da states (integrated only over the electron angular 
coordinates) behaves as at r — > 00. To account such 
a behavior in the framework of ordinary (Coulomb field 
only) MQDT, one would have to use energy-dependent 
non-diagonal matrix element of the quantum defect ma- 
trix. For this reason, we cannot use the specified basis of 
electronic states to represent the electronic Hamiltonian 
of the system. However, we adopt the logic of f23l. [2^. [25| 
and utilize a different electronic basis for which the non- 
diagonal coupling elements would vanish at long range 
(or at least decrease with r faster than I/r^). 

We stress the two approximations used here: (1) We 
consider the dipole moment for all possible geometries of 
HCO''" ion to be constant and equal to the permanent 
dipole moment D of HCO"*" at its equilibrium position 
(linear configuration, 9 = 0). (2) As we already men- 
tioned, da states have a minimal influence on the DR 
process and, therefore, we can neglect any coupling with 
the states and concentrate our attention uniquely on the 
coupling between pa and sa states. In fact, we have ver- 
ified the validity of the second approximation by com- 
paring results with and without inclusion of couplings 
between pa and da states, and we find that the inclu- 
sion of the coupling with the da states has a negligible 
effect on the final cross-section. To enhance the clarity 
of our presentation, we do not include the coupling with 
the da states in the discussion below. With the men- 
tioned approximations, the sa-pa part of the electronic 
Hamiltonian Hei{Q,r) has the form 

/ _1 <Y'|°|£)cos6)e|y°> \ 



(1) 

The above matrix is coupled to the two pn^ states by the 
Renner-Teller coupling between pir and pa states. The 
matrix elements for the coupling are given and discussed 
inRefs. 

Under the assumption that D is constant for any Q, 
the Hamiltonian depends only on r. The numerator in 
the non-diagonal elements is easily evaluated and we find 
D = {Y^\Dcosee\Y^) = ^ (here the angle 6^ is the az- 
imuthal angle of the electron in the molecular coordinate 
system). When we diagonalize the 2x2 Hamiltonian 
above, we find two new electronic states \sa) and \pa), 
which are superpositions of \sa) and \pa) with a projec- 
tion of the angular momentum still equal to zero. The 
radial functions related to these new channel eigenstates 
are solutions of two Schrodinger equations with effective 
potentials different than for \sa) and \pa) states: 
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V±{r) 



1± V1 + 4D2 



2r2 



(2) 



The potentials V+ (r) and VL (r) are respectively related 
to the channels \pa) and \sa). The effective potential for 
\'pa) is more repulsive than the one for \pa) ; the effective 
potential for \sa) is more attractive than the one for \sa). 
We use the notations \pa) and \sa) for the new channel 
states because the mixing between pu and sa states is 
not very big: the new states are still mainly composed 
from their respective original states (around 75%). 

Changing electronic states from \pa), \s(j) to \p(j) and 
\sa), the 5x5 electronic Hamiltonian H = Hei{Q) in- 
cluding the five states \pa), \sa), \da), and \ptt^) has the 
same form as the one given in Refs. Q if we change 
quantity on Refs. [s, f6l| to E^ to keep consistency in 
notations. The coefficients 7 and 6 are obtained from the 
ab initio potential energy surfaces VTr\Tr",a{Q) using 
the same formulas (3), (5), and (6) in Ref. where the 
states denoted as pir' and ptt" , are respectively symmet- 
ric and antisymmetric superposition of Iptt"*") and |p7r~) 
with respect to reflexion in a plane containing the molec- 
ular axis. Therefore, the diagonalization oi H = Hei{Q) 
is accomplished by the same unitary matrix U as in Eq. 
(2) of Ref. g. 



III. QUANTUM DEFECTS AND REACTION 
MATRIX FOR THE -l/r±A/r'^ POTENTIAL 



the following equation, where A is a complex number: 
1 (P A(A + 1) 1 



2dr2 



2r2 



e /(e,A,r)=0. (4) 



The equation has two independent solutions behaving 
asymptotically for e < as /*(e, A,r) e^^r^'^, where 
1/ = ^ = (— 2e)~^/2 (Pqj. j± definitions of Ref. 

[3. See Eq. (3.7) of [l6|]. In contrast, /± in Ref. [H are 
defined differently, Eq. (2.38) in [13,]). From these two 
functions the regular and irregular solutions for this equa- 
tion (respectively denoted /(e, A, r) and g{e, A, r)) can be 
constructed. In our case, the only basic requirement for 
g{e, A, r) is that it should lag 90° with respect to the reg- 
ular solution /(e. A, r). More generally, great care has to 
be taken to construct the irregular solution for certain 
types of potentials (as pure dipole or repulsive Coulomb 
potential) in a way that the solution remains indepen- 
dent on threshold effects. This construction procedure is 
discussed for instance in Ref. [lB|- The asymptotic ex- 
pansion of these solutions as a superposition of is ob- 
tained using the following general formula. For instance, 
for /(e. A, r) we have: 

2zfc/(e,A,r) = 
W^(r, /)/+(£, A, r) - Wif+, f)f-{e, A, r) , (5) 

where the Wronskian W{a, b) above is defined as 
W{a,b) = - b^. For the potential in Eq.®, the 
solutions are confluent hypergeometric functions. Using 
their asymptotic behavior at large r, we can write for / 



One more difference with the treatment of Refs. [1, Q 
is the way how the quantum defects functions fJ.i{Q) and 
the reaction matrix K{Q) are evaluated from the ab ini- 
tio potential surfaces of HCO. The evaluation of quan- 
tum defects ^i(Q) in Refs. d, 0, [2^ is made based on 
the Rydberg formula, which assumes an integer and non- 
negative partial wave quantum number and the asymp- 
totic effective potential for the electron-ion interaction of 
the form V{r) — + ^'"2^2^'^ with A an integer usually 
denoted/. From Eq. ([2]) it is clear that the corresponding 
values of A are not integer numbers: they are obtained 
from the equation A(A +1) = 1 ± \/l + 4JD^ . Intro- 
ducing the positive A_|_ = 5 + 4vl + 4~D2 and negative 

A_ = h — Ayl + AEP quantities, the new "partial wave" 
quantum numbers are 



Apcr 



_ -l±i^-A_ 
2 



(3) 



Apcr is real, Xsa is complex. Therefore, the standard quan- 
tum defect treatment cannot be applied. Below, we sum- 
marize the main ingredients of the generalized quantum 
defect theory as formulated in [l^ [3] ■ 

We start with an one-channel problem represented by 



/(e,A,r) = {-K n)-^- {sin (3{D,)-\r -cos I3DJ+). (6) 

In the above expression, the factor is a monotonic 
function of e that rescales the amplitudes of and 
the oscillating factors sin (3 and cos (3 are parameters re- 
sponsible for the energy-dependent mixing between the 
exponentially growing and decaying functions The 
wave function for a bound state should decay exponen- 
tially at infinity. Therefore, the condition for a bound 
state is simply sin/3=0 or equivalently (3—TT(n + 1) with 
n = 0, 1.2.... If A is a real number, /3 is simply 7r(i/ — A). 
It corresponds to a pure Coulomb field with a centrifugal 
term and we recover the Rydberg formula for the energy 
of bound states E = —0.5/ {n + l + 1)^. In a general case 
when A = Xr -I- ia is a complex number, the formula for 
/3 is more cumbersome [l6j 

/?(«:, A) = 

7r(l/K — Xb.) + tan""'" [tanh(7rQ;) tan(y — a ln(2K)] , (7) 
where y = arg(r(2A + 2)[r(i^ - A)/r(A + 1 + ly)]^/^) . 

Having defined the regular solution, we can write the ir- 
regular solution using a similar procedure. The irregular 
solution g{e, A, r) is then given by 

gie,X,r) ^ -i7rK)-^cospiD,)-'r +sinl3DJ+) . (8) 
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If we consider an additional short range potential in 
Eq. (HI, the asymptotic solution will in general be writ- 
ten as a superposition of regular and irregular solutions 
F(e, A, r) — /(e. A, r) — Kg{e, A, r), where if is a mixing 
coefficient. From Eqs. ([6]) and ([8]) the asymptotic behav- 
ior of F{e,X,r) is easily derived. The requirement that 
the coefficient in front of the growing exponential must 
vanish in order for the solution to be physically accept- 
able determines K as K = — tan(/3). 

Because £[00+ has a dipole moment, evaluation of 
the reaction matrix K from ab initio energies has to be 
changed (in comparison with our previous study 
Namely, the reaction matrix K in the diagonal form and 
at fixed nuclear positions is written in terms of /9(k, A) as 
Kii = — tan{Pi{Ki, Xi)) for each channel i (we have five 
of them as discussed above). Index i at Pi means that 
we use different formulas depending on the asymptotic 
behavior of the potential in the corresponding channel. 
The value of k is obtained from the ab initio electronic 
energies of HCO and HCO"'" [2^ with n ~ A for sa state 
and n = 3 for other states. The K{Q) matrix in the basis 
of states StJ, pvr", per, pTr"*" and da is obtained using the 
unitary transformation as discussed in Ref. [1]: K{Q) = 
J/t tan(-/3)C/. 



+ experiment by Le Padellec et at. [28] 
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FIG. 1: (Color online) Calculated DR cross section for HCO^ 
including dipole contribution (solid line) and the DR cross 
section for HCO'*' without dipole contribution as a function 
of the incident electron energy. The experimental p7j (cross 
symbols) and previous theoretical (dashed line) cross sec- 
tions are also shown for comparison. The theoretical curves 
include a convolution over the experimental electron energy 
distribution. 



IV. CROSS-SECTION CALCULATION 

Similarly to our previous study, we solve the vibra- 
tional Schrodinger equation for IICO+ fixing the adi- 
abatic coordinate i?GH and obtain the eigenenergies 
Um^,i{RGn) (we call them adiabatic potentials) and the 
corresponding eigenstates $m^,i(i?GH; .Rco, ^, <p) (adia- 
batic states). Then, we use $m^,i(i?GH; RcOj^: f) to ob- 
tain the reaction matrix /C{„^,i^i},{,„/^,i/,i/} (-Rgh) 

(9) 

where the integral in the above expression is taken over 
the three coordinates, i?cOj and 9. The resulting re- 
action matrix ICjj' (j = {mip,l,i}) is multichannel and 
depends parametrically on Ron- It is then used to ob- 
tain for each Rqh value positions UaiRcii) and widths 
^aiRau) of autoionizing resonances of the neutral HCO 
molecule. The procedure of resonance search is standard 
and is discussed, for example, in our earlier study The 
only difference is how the quantity (3{k, A) is calculated 
for states with non-integer A. (3{k, A) is needed to per- 
form the closed-channel elimination procedure and eval- 
uation of closed-channel mixing coefficients. (For more 
details about the channel elimination and the closed- 
channel mixing coefficients see Eqs. (2.50) and (2.54b) 
and the corresponding discussion in Ref. [13.]). The 
quantity (3{k, A) should be calculated using Eq. ([7]) for 
channels with non-integer A. The obtained widths and 
positions of autoionizing resonances are used to calculate 
the electron-ion recombination cross-section as described 
in Ref. 5]. 



V. DISCUSSION OF RESULTS 

Figure[l]shows the cross section obtained in the present 
study. The inclusion of the electric dipole moment of 
HCO'*" leads to an overall increase of about 50% for the 
cross section and agrees better with experiment than the 
cross section calculated neglecting the dipole moment. 
The agreement between theory and experiment is very 
good for electron energies between 0.06 eV and 0.3 eV. 
For larger energies, the experiment seems to have a large 
error bar (the experimental error bar is large possibly 
because the measured cross-section is small at higher en- 
ergies). At energies below 0.06 eV the theoretical cross- 
section is below the experimental one by a factor of two 
approximately. The reason for the difference is not clear. 
One possible reason for the remaining disagreement could 
be related to the isomer HOC+ of the HCO+ if (1) the 
isomer is present in the experiment of Ref. [27j in a non- 
negligible fraction and (2) the DR cross-section for the 
isomer is significantly larger than the cross section for 
HCO"*". The preliminary calculation [28,] suggests that 
the ground ionic potential of HOC+ crosses the dissocia- 
tive potential of HOC in the Franck-Condon region. It 
means that the DR cross section for HOC+ should in- 
deed be larger and could be of order of a few 10~ ''cm^/s 
at 0.025 eV (^300 K). Nevertheless, this hypothesis must 
be verified in experiment and calculation. 

We would like to discuss in some detail the increase of 
the cross section comparingwith the results of the pre- 
vious theoretical approach Q . With the effect of the 
HCO"*" electric dipole moment included, it seems reason- 
able to expect a higher cross section due to an additional 
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coupling between sa and pa states introduced by the 
dipole term. However, a more detailed analysis shows 
that the coupling is not directly responsible for the in- 
crease of the cross section, but it acts rather in an indirect 
way through the Renner- Teller effect as described below. 
The cross section in the present study is calculated us- 
ing autoionizing resonances. Therefore, the cross section 
is increased if the density or/and widths of resonances 
become larger with the inclusion of the electric dipole 
moment. The density of resonances is expressed as ^|§. 
Therefore, we can compare the density of resonances by 
evaluating this quantity using the updated formula for f3. 
For the p-states, the quantities f3 for the pure Coulomb 
and Coulomb-|-dipole interactions differ from each other 
by an energy-independent term. Therefore, the density 
of resonances produced by closed p-wave channels is the 
same in the both treatments. For the s-states, analysis 
of the derivative of new /3{E) of Eq. ([7]) gives only negli- 
gible differences for pure Coulomb and Coulomb-f dipole 
interactions. Thus, the higher observed cross section is 
caused uniquely by an increase in widths of autoioniz- 
ing resonances. By inspecting different resonances, we 
have indeed found that the increase in the cross section 
is related to larger widths of p- resonances. To apprehend 
this result, we consider non-diagonal elements of the K 
matrix. These terms can be written in the following way 

K^AQ^E)=-7r<q^.,E\V\4,,>E> + 

w I dE^^tmmpE:^K,,[Q^E'), (10) 

where indexes «, i', and k specify channels with the elec- 
tronic wave function 4>, the symbol V refers to the prin- 
cipal value of the integral, and V is the potential of the 
system. Now in a first approximation, we can just retain 
the first term on the right of Eq. As the width of 

a Feshbach resonance is 

dE 

V^2tt—\<4>:e\V\4>.,'e>\\ (11) 
where n is the principal quantum number, we obtain 

T^—^\K^.{Q,E)\'' . (12) 

Because the widths are related to the non-diagonal el- 
ements of the reaction matrix K^-i {Q, E), we can inves- 
tigate the non-diagonal elements to understand the rea- 
son for the increase of the cross section. In the diag- 
onalized form, the elements of the reaction matrix are 
Kjji = tan(7r/ij)(5jj' where the quantum defects fj,j are 
obtained directly from ab initio calculation for every in- 
ternuclear configuration Q using Eq. ( [7]) . The old quan- 
tum defect are also obtained from the same ab initio data 
but using the Rydberg formula. The new and old quan- 
tum defect are shown for comparison in Fig. [5] as a func- 
tion of the Jacobi angle 6. As it can be seen from the 
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FIG. 2: (Color online) The figure shows quantum defects for 
the new states sa and pa calculated from the ab initio elec- 
tronic energies using l3 of Eq. ((Tjl as well as the old quantum 
defect (calculated using /3 determined by a pure Coulomb po- 
tential) as a function of the Jacobi angle 6. The two other 
Jacobi coordinates -Rco ~ 2.00 a.u. and Ran = 3.27 a.u are 
fixed at their equilibrium positions for 9 — 0. 



figure, the new quantum defect for the sa and pa states 
differ from each other by a simple translation. The trans- 
lation, which has no direct effect on widths of the sa 
resonances (they are basically just shifted), has a drastic 
effect on the pa states coupled to the pir states by the 
Renner-Teller effect. It is possible to show that changing 
the basis for the reaction matrix from the eigenbasis (di- 
agonal form of K{Q)) to the basis of the five states \pa), 
\sa), \da), and |p7r*) discussed above, introduces a non- 
diagonal term Kui between pa and pn^ that is directly 
proportional to the difference tan(/ipCT) "tan(/ip^' ). From 
Fig. [21 it is evident that the difference is increased for 
the new pa state comparing with the difference for the 
old (without dipole contribution) pa state. The coupling 
element between pTr'*' and p'K~ are only slightly affected 
by the translation of the quantum defect. 

The approximate width for the two coupling elements 
in the new and old can be estimated using Eq. . The 
result is shown in Fig. [3l The coupling between pa and 
pir"^ is clearly increased compared to the previous treat- 
ment by about 50%. On the other hand, the coupling 
between the p7r+ and p'K~ states is changed only a little 
(less than 5%). Thus, we conclude that the increase of 
the coupling between pa and pn^ states is responsible for 
the higher cross-section in the present approach includ- 
ing the electric dipole moment of IICO+. We see also 
that such a rough estimation of the increase in resonance 
widths is compatible with the overall increase in the fi- 
nal cross section obtained in the full calculation, which 
proves the validity of the interpretation. 
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tion of the ion with low energy electrons. We found that 
the inclusion of the dipole moment increases significantly 
the coupling between the pa and pn^ states and only a 
little the coupling between pTr"*" and p7r~. This leads to 
an increase of the cross section for all collision energies 
below 1 eV by about 50%. The obtained theoretical cross 
section agrees well with the experimental data of Ref. 
[2?! for energies between 0.06 and 0.3 eV. Although, the 
cross section is also increased for low energies (< 0.06 
eV), the theoretical result is still below the experiment 
by about a factor of two. More theoretical and experi- 
mental work has to be done to clarify the reason for the 
remaining difference. 



FIG. 3: (Color online) The figure shows the absolute value 
square of the non-diagonal elements \K^^i{Q,E)'^ of the re- 
action matrix as a function of the bending angle 6 around 
equilibrium values of 7?co(2.00 a.u) and Ran (3.27 a.u). 



VI. CONCLUDING REMARKS 

We have investigated the role of the permanent electric 
dipole moment of HCO"*" in the dissociative recombina- 
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